Kinetic energy sum spectra in nonmesonic weak decay of 

hypernuclei 

Cesar Barbero^''^, Alfredo P. Galeae^, Mahir S. Hussein^'^, and Francisco Krmpotic^'^'' 
^ Facultad de Ciencias Exactas, Departamento de Fisica, 
Universidad Nacional de La Plata, 1900 La Plata, Argentina 
^ Institute de Fisica Teorica, Universidade Estadual Paulista, 
Rua Pamplona 145, 01405-900 Sao Paulo, SP, Brazil 
^Departamento de Fisica Matemdtica, 
Instituto de Fisica da Universidade de Sao Paulo, 
Caixa Postal 66318, 05315-970 Sao Paulo, SP, Brazil 
^Instituto de Fisica La Plata, CONICET, 1900 La Plata, Argentina 
^ Max-Planck-Institut fiir Physik komplexer Systeme, 
Nothnitzer Strafie 38, D-01187 Dresden, Germany and 
^Facultad de Ciencias Astronomicas y Geofisicas, 
Universidad Nacional de La Plata, 1900 La Plata, Argentina 
(Dated: October 17, 2008) 



1 



Abstract 

We evaluate the coincidence spectra in the nonmesonic weak decay (NMWD) A7V — nN of A 
hypernuclei ^He, aHc, j^^C, O, and ^^Si, as a function of the sum of kinetic energies EnN = 
En + Ejsi for A'^ = n,p. The strangeness-changing transition potential is described by the one- 
meson-exchange model, with commonly used parameterization. Two versions of the Independent- 
Particle-Shell-Model (IPSM) are employed to account for the nuclear structure of the final residual 
nuclei. They are: (a) IPSM-a, where no correlation, except for the Pauli principle, is taken into 
account, and (b) IPSM-b, where the highly excited hole states are considered to be quasi-stationary 
and are described by Breit-Wigner distributions, whose widths are estimated from the experimental 
data. All np and nn spectra exhibit a series of peaks in the energy interval 110 MeV < E^n < 170 
MeV, one for each occupied shell-model state. Within the IPSM-a, and because of the recoil effect, 
each peak covers an energy interval proportional to A , going from ^ 4 MeV for f Si to ^ 40 MeV 
for jjie. Such a description could be pretty fair for the light ^He and ^He hypernuclei. For the 
remaining, heavier, hypernuclei it is very important, however, to consider as well the spreading in 
strength of the deep-hole states, and bring into play the IPSM-b approach. Notwithstanding the 
nuclear model that is employed the results depend only very weakly on the details of the dynamics 
involved in the decay process proper. We propose that the IPSM is the appropriate lowest-order 
approximation for the theoretical calculations of the of kinetic energy sum spectra in the NMWD. 
It is in comparison to this picture that one should appraise the effects of the final state interactions 
and of the two-nucleon-induced decay mode. 

PACS numbers: 21.80.-|-a, 13.75.Ev, 21.60.Cs, 21.10.Pc 

Keywords: nonmesonic hypernuclear decay; energy spectra; strength functions; one meson exchange model 
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I. INTRODUCTION 



Because of the difficulty in detecting neutrons, up to a few years ago only the high 
energy proton spectra had been measured in the two-body AA^ —>■ nN nonmesonic weak 
decay (NMWD) of A hypernuclei. The proton-induced transition rates Tp = r{Ap np) in 
|He, \^Li, and have been determined in this way . The corresponding neutron- 

induced transition rates r„ = r{An —>■ nn) were estimated through the comparison of the 
measured proton spectrum with that of the intranuclear cascade (INC) calculation where the 
r^/Fp ratio is treated as free parameter. This procedure signalized very large experimental 



ireeparE 

m 



n/p ratios (= 0.9-2.0) [l|, 1^, l3|, in comparison with theoretical estimates 



10 



0.3 - 0.5) 



nj. (See, 



obtained within the One-Meson-Exchange Potential (OMEP) 

for instance, 0, Table 3].). In spite of large uncertainties involved in the indirect evaluation 
of neutrons, this discrepancy between data and theory was considered to be a serious puzzle 
in the NMWD. 

Yet, quite recently, the above scenario has drastically changed due to the very significant 
advances in our knowledge (mainly experimental) on both neutron and proton spectra. 
These advances are: 

1. In the experiment E369 were measured high-quality neutron spectra in the decays of 
^^C and ^^Y, which made possible to compare them directly with the corresponding 



A ^ '^'■'^^ A 

proton spectra, yielding the result r^/Fp = (0.45 — 0.51) ± 0.15 



121. 



2. Okada et al. [13] have simultaneously measured the energy spectra of neutrons and 
protons in ^He and at a high energy threshold (60 MeV), from where it was 
inferred that for both hypernuclei F„/Fp = 0.5. 



3. Garbarino, Parreno and Ramos [l6j have called attention on the fact that 



lan single 



13, 
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"correlation observables permit a cleaner extraction of F„/Fp from data t 
nucleon observables", which has stimulated several experimental searches 
I20I ]. They have also done a theoretical evaluation of the pair distributions for: a) the 
sum of the nN kinetic energies E = En + En, SnN^E), and b) the opening angle 9, 
SnN{cos9), in the decays of \B.e and j^C. The number of detected nN pairs N„Ar is 
proportional to F^v = / SnN{E)dE = J SnN{cos9)dcos9, and therefore N„„/N„p = 
F„/Fp (see below). The primary weak decay was described by the OMEP dynamics 



within the framework of a shell model, while an INC code was used to take into account 
the strong final state interactions (FSI) involving the primary nucleons and those in 
the residual nucleus, including the possibility of emission of secondary particles. They 
conclude that the datum N„„/N„p = 0.44 ± 0.11 obtained in KEK-E462 ^| for |He 
is compatible with an n/p ratio of 0.39 ±0.11 for this hypernucleus if the two-nucleon 
induced mode ANN nNN is neglected, or an even lower value if it is included. 



2l|, they found Tn/TpQ^C) = 0.46 ± 0.09 when the 



Similarly, together with Bauer 
two-nucleon induced mode is neglected and a slightly lower value when it is included. 
The shell-model used in Refs. jl4, 1^ is substituted in the last work by a nuclear 
matter formalism extended to finite nuclei via the local density approximation, while 
the Monte Carlo INC model is retained to account for the FSI. 



4. Quite recently Bauer [22] has given a step forward with his nuclear matter formalism, 
describing microscopically both the weak decay mechanism and the FSI, confirming 
in this way that the latter lead to somewhat lower value for the ratio F^/Fp in ]^C. 

5. Several coincidence emission measurements of the above mentioned spectra have been 
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20|, 



23 



2J], from which were extracted the new experimental 



performed p/7|, ll 

resuhs, F„/Fp(iHe) = 0.45 ± 0.11 ± 0.03 and F„/Fp(jv^C) = 0.40 ± 0.09, that point 
towards the solution of the longstanding F„/Fp puzzle.^ 



Reasoning within the two-body kinematics for the one-nucleon induced NMWD, it is 
expected that SnN{E) should exhibit a narrow peak at the two-particle energy E close to the 
decay Q- value, while SnN {cos 6) should be restrained within the back-to-back angle 6 = it. 
Thus, all experimentally observed deviations from such spectral shapes are very frequently 



3,y 



23 



attributed to the FSI and/or to the two-nucleon induced processes |l7l . 
However, this is no longer the case when the recoil of the residual nucleus is taken into 
account, which makes the kinematics to be of a three-body type. Moreover, when the shell 
model structure is also taken into account the energy spectra will have a bump at each 



^ We note that the relationship r„/rp(^He) ~ Tn/^p(j^C) could have a very simple explanation, similar 
to that given iri Ref. [25|] for the asymmetry parameter: aA(^He) ~ a\(j^G). 



2 In Refs. 



ven iri Kei 
M lli it i 



is said that the np energy spectra of the three-body proton-induced decay ^He H 



+n + p should exhibit a narrow peak close to its Q- value of 153 MeV, which is only valid when the recoil 
effect is neglected. 
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single particle state, the width of which will depend on the magnitude of both the recoil 
and the spreading in strength of the hole states in the inner shells. Although the detailed 
structure and fragmentation of hole states are still not well known, the exclusive knockout 
reactions provide a wealth of information on the structure of single-nucleon states of nuclei. 
Excitation energies and widths of proton-hole states were systematically measured with 
quasifree (p, 2p) and (e, e'p) reactions, which revealed the existence of inner orbital shells in 



nuclei 
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21, 



28 



29 



3G, 



31 



32 



33 



34|. 



The aim of the present work is to discuss quantitatively the interplay between the recoil 
effect and the nuclear shell structure in the kinetic energy sum spectra of NMWD. The paper 
is organized as follows. In Section [III we discuss the calculation of these spectra within the 
Independent-Particle-Shell-Model (IPSM). In Section [XTTl we exhibit the numerical results for 
^He, iHe, ]^C, jfO, and fSi hypernuclei. In Section [IV] we discuss these results and their 
connection with the experimental data. Finally, in Section |Vl we present several concluding 
remarks. 



II. KINETIC ENERGY SUM SPECTRUM 

For the purpose of completeness and clarity, we will first redo the calculation of the decay 
rate Fjv using Fermi's golden rule. The novelties here, in comparison with our previous works 
Refs. p, 0, [sj, are the recoil effect and the spreading of the single-particle configurations. 
As will be shown in Section IIV| their role is of minor importance in the evaluation of the 
integrated transition rates Fjy, as well as on the ratio F„/Fp, but they are crucial for a correct 
description of the energy distribution of the transition strength. In fact, a single-particle 
state Ijat) that is deeply bound in the hypernucleus, after the NMWD can become a highly 
excited hole-state IJ^^) in the continuum of the residual nucleus. There it suddenly mixes 
up with more complicated configurations (2hlp, 3h2p, . . . excitations, collective states, etc.) 
spreading its strength in a relatively wide energy interval {sS^.^ This happens, for instance, 
with the lsi/2 orbital in ^^C, that is separated from the lp3/2 state by approximately 23 
MeV, which is enough to break the 10 particle system, where the energy of the last excited 
state amounts to ~ 16.5 MeV. 

^ One should keep in mind that the mean hfe a A hyperon is ta = 2.63 x 10^^" s, while the strong interaction 
times are of the order of 10^^^ s. 
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The NMWD rate of a hypernucleus (in its ground state with spin J/ and energy Ej^) to 
residual nuclei (in the several allowed states with spins Jp and energies Ea^jp) and two free 
nucleons nN (with total spin S and total kinetic energy E^n = En + E^), reads 



27r J\{PnPNSMs;ar,JFMF\V\JjMi) 



N 

SMsonJfMp 

X d{A^^jp-Er-EnN)j^J^. (1) 

Here, V is the hypernuclear nonmesonic weak transition potential, and the wave functions 
for the kets IpuPnSMs; unJeMf) and jJ/M/) are assumed to be antisymmetrized and 
normalized. The label stands for different final states with the same spin Jp, E^ is the 
recoil energy of the residual nucleus, and 

i^^^jp = A + Ej^-E^^jp, with A = Ma-M = 176 MeV, (2) 

is the liberated energy. The two emitted nucleons are described by plane waves, and initial 
and final short range correlations are included phenomenologically at a simple Jastrow-like 
level. 

It is convenient to perform a transformation to the relative and cm. momenta (p = 
\{Pn — Pn), P = Pn + Piv), coordinates (r = r„ — r^v, R = |(r„ + riv)) and orbital angular 
momenta 1 and L, and to express the energy conservation as 

EnN + Er — A^jvJf = + ~ ^a^Jp = 0, (3) 

where 

'^~M' '~ 2M{A-2y '''~AM7^~2 '' ^ ' 

are, respectively, the energies of the relative motion of the outgoing pair, of the recoil, 

and of the total cm. motion (including the recoil). Following step by step the analytical 



developments done in Ref. 



6|, the transition rate can be can expressed as 



N 



lo dtp 

where we have defined (see [S, Eqs. (2.13) and (2.14)]) 



dep (5) 



'^^P Tl \ A / SXILTJcnJf 



X 



J2M{pPILXSJT;jn){Ji\\ {a}^cil) Jla^Jp) 



JN 



2 



(6) 



6 



it being understood that the square root is to be replaced by zero whenever its argument is 
negative. The angular momentum couplings I + L = A and A + S = J have been carried 
out, J = \J2J + 1, and y4 = Z + A^ + lis the total number of baryons. 

It is self-evident that for A oo one obtains the same result as in Refs. . It is 

also worth noting that the overall outcome of the recoil on Fat is very small, mostly because 
the effect of the factor (^^) ^ in Eq- dS]) is, to a great extent, cancelled by the effect of 
the factor (^442) originating from ^Jep{K^^^~j~-7p)dep. This is the reason why we have 
not included the recoil previously. 

The spectrum of Fat as a function of EnN is now easily obtained from Eq. ([6]) by means 
of the relation 

2 

Calling E = EnN, one gets 

A 





as follows from ([3]) and 



N 



SnwiE) dE 



(8) 



with 



SnN{E) 



TT 



X 



SXlLTJajsi Jjr 

2 



Y,M{pPlL\SJT-3M){Ji\\ {a]^a]^)j\\aMJF) 

3n 



(9) 



where 



P 
P 



I MA 



E-A' 



afjjp j 5 



A' 



2M(A-2)(A„^j^-E), 
A-2 



Q-nJf 



A 



and the condition 



A' 



cinJf 



<E < A 



A 



(10) 

(11) 

(12) 



has to be fulfilled for each contribution. 

As previously , it will be assumed that the hyperon in the state j\, with single- 

particle energy e^^, is weakly coupled to the A — 1 core, with spin Jq and energy Ec = 
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Ejj — ej^. Then the initial state is | J/) = \{JcjA)Ji) : and the spectroscopic amphtude 
(J/|| (ctjjv'^L) J W'^nJf) can be rewritten as 

(JiWialaD^Wa^Jp) = JJj T"" \ {Jc\\al\\a^Jp). (13) 

[ J 3n Jf J 

The following two approaches for the final states \aNJF) will be examined within the 
IPSM. 



A. IPSM-a 

Here, we completely ignore the residual interaction and, consequently, the only states 
IctArJi?) giving a nonzero result in Eq. f|T3|) and therefore contributing to Eq. ([9]) are those 
obtained by the weak coupling, and properly antisymmetrizing, of the one hole (Ih) states 
|j^i) to the core ground-state \Jc)- That is, recalling that we completely ignore the residual 
interaction in this approximation, 

WnJf) ^ \3nJf) = \{Jc,3n'^)Jf) , and K^^J^. ^ = Eq - ej^ , (14) 

where e^^ is the single-particle energy of state j^- As an illustration, in the case of |^Si the 
model space contains four single-particle states, both for protons and for neutrons (n^ = 
n„ = 4), namely, lsi/2, 1^3/2, lpi/2 and 14/2- Thus, for \ Jc) = 114/2'^-""'^), the final states 
(fT4|) are constructed by adding two holes in the "^^Si nucleus, and read: 

f Si ^nn + 26Si fSi ^np + ^^A\ 



(15) 



1(14/2^-1)2; 0, 2, 4) 1(14/2^-114/2^-'); 0, 1, 2, 3, 4, 5) 

|14/2ra-ilsi/2ra-i; 2, 3) |14/2ra-ilsi/2p-i; 2, 3) 

|14/2n"il]9i/2n"i; 2, 3) |14/2n-ilpi/2]9"i; 2, 3) 

|14/2ri-ilp3/2ri-i; 1, 2, 3, 4) 1 14/2^^-1 1P3/2P-^ 1, 2, 3, 4). 

The summation on Jp in iQ can be performed for each single-particle state Jn, as done 
in 0, Eqs. (11), (12), (13)]. One gets 

4M^ I ; 

SnN{E) = -—^A{A - 2)3 5: ^(A,, - E){E - A^-J^.^pP), (16) 
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where 

J=iiV+l/2 

^UpP) = FjiT. M\pPlLXSJT;j^), (17) 

J=|jjv-1/2| SILXT 

with the spectroscopic factors exhibited in 0, Table 1]. The maximum and minimum 
hberated energies are, respectively, 

A,, J, ^ A,^ = A + e,^ + e,, (18) 

and 

= (19) 

and the momenta are given by 
P = 



B. IPSM-b 



I MA 



(E-Ag and P = ^2M{A-2){A,,-E). (20) 



Formally, one starts from the unperturbed basis I^at Jf)o with i^v = 1, 2, . . . n^, n7v + 1, • • •, 
where for za? < we have the same simple doorway states IJnJf) in Eq. f|T^ (listed in 
Eq. (ITS!) for |^Si), while for i^v > n^v + 1 we have more complicated bound configurations 
(such as Ship, 4/i2p, ... in the case of A^Si) as well as those including unbound single-particle 



states in the continuum. As in Ref. [35|], the perturbed eigenkets IonJe) and eigenvalues 
Ea^jp are obtained by diagonalizing the matrix Q{i]<qJp\H\i'j^Jp)Q of the exact Hamiltonian 
H: 

(21) 

with 

oo 
iN=i 

oo 

= T.C-:''njNJF)+ E Cr:'n^NJF)0. (22) 

jjv i]V=njv+l 

It is easy to see that only the ket IJnJf) in the expansion fl22l) will contribute to the matrix 
element (Jc||a]^||aArJF) in Eq. (fT3|) . Therefore Eq. (Q takes the form 



71 * 

■ya jv Jf I 2 



X E IV(^"-^- - ^)(^ - K^J.)^(p^P-^JnJf) , (23) 



where 

Hv.P-.JnJf) = JT^ E \M{pPlL\SJT-3M){Ji\\{a]^a]^)j\\3NJF)^ . (24) 



ILXSJT 



To evaluate the amphtudes C*"^"'^ one would have to choose the appropriate Hamiltonian 
H and the unperturbed basis \iNJF)o, and solve the eigenvalue problem (1^ . We will not 
do this here. Instead, we will make a phenomenological estimate. First, because of the high 
density of states, we will convert the discrete energies Aq,^j^ into the continuous variable e, 
and the discrete sum on a^r into an integral on e, i.e., 

roo 

A«,Jf ^e, E \Cl-'-\' ^ E / \C,,j,ie)\'pjAe)de, (25) 

a^jp Jp 

where pjp{e) is the density of perturbed states with angular momentum Jp. In this way the 
spectrum ( l23l) becomes 

4M3 



SnNiE) = V^(A-2)3 

/OO I 
P..vJf(^)V -E){E- e')J^{p, P;jNJF)de , (26) 
-OO 



X 



where 

P,^jpie) = \C,,j,ie)\'pj,{e) (27) 



is called the strength function 35|, |36|, |37|] and represents the probability of finding the 



configuration IJnJf) = \{.Jc^ 3n^)Jf) per unit energy interval. Moreover, 



and the condition 



P = ^2M{A-2){e - E), (28) 
e' = e^-, (29) 



e' <E <e, (30) 



has to be fulfilled throughout the e integration. It is convenient to introduce the averaged 
strength function 



Jc+Jn 

E 

JF = \Jc-jN\ 
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where 



dim(jArJc) 



2jN + 1 for jjv < J( 



c 



2Jc + 1 for Jc < JN ■ 

This aUows to simphfy Eq. fl26l) by making the approximation Pjj^jp{£) ~ Piiv(^) to get 



(32) 



4M^ I 

SME) = a/A(A-2)3 

TT * 



5: r P,,ie)^ie-E)iE-e')J^j,ipP)de 

J —oo 

IN 



(33) 



where we noticed that the summation of J^{p, P; JnJf) over Jp gives the function defined 
in Eq. (fTT]). 

The IPSM-a results would be recovered if one made the further approximation 



P,,(e) = 5(£-A, 



JnJ- 



(34) 



Here, in IPSM-b, the 5-functions ( IMI) will be used for the strictly stationary states, while 
for the fragmented hole states we will use Breit-Wigner distributions, 

27,^. 1 



71 7? +4(£-A, 



(35) 



where 'jj^ are the widths of the resonance centroids at energies Aj^ (see [35!, Eq.(2.11.22)]). 

It might be important to point out that, since both strength functions Pj^{e) are nor- 
malized to unit, their effect on integrated observables like the decay rates Pat is expected 
to be small even if they considerably affect the spectra. This will be further investigated in 
Section IVl 



III. NUMERICAL RESULTS 



In Figs. [U and [2] we show, respectively, the normalized energy spectra Snp{E)/Tp and 
Snn{E)/Tn for ^He, iHe, ]^C, fO, and i^Si hypernuclei, evaluated within the full OMEP, 
that comprises the [71,1], K, p,uj, K*) mesons. The single-particle energies for the strictly 



stationary hole states have been taken from Wapstra and Gove's compilation 



38| , and those 



of the quasi-stationary ones have been estimated from the studies of the quasi-free scattering 



processes {p, 2p) and (e, e'p) 



26 



27 



28 



29 



30 



31 



32 



33 



34| 



The two IPSM approaches exhibit some quite important differences: 
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FIG. 1: (Color online) Normalized energy spectra Snp{E) /Tp for ^He, ^He, a^C, \^0, and ^*Si hypernuclei 
for the full OMEP obtained within the approaches IPSM-a (upper panel) and IPSM-b (lower panel). For 
the s-shell hypernuclei, only the IPSM-a approach has been used. 

a) IPSM-a: The spectra cover the energy region 110 MeV < E < 170 MeV and contain 
one or more peaks, the number of which is equal to the number of shell-model orbitals 
lsi/2, lp3/25 lpi/25 1*^5/25 2si/2, 1^3/2 ' ' ' that are either fully or partly occupied in \ Jc)- 
Before including the recoil, all these peaks would be just spikes at the liberated energies 
Ajj^, as can be seen from ^ setting Er = 0. With the recoil effect, they behave as 

SME = A,J ~ ^(A,, - E){E - A;.Je-^■^(^-2)(^--^)^^ (36) 

and develop rather narrow widths ~ [b'^M{A — 2)]^^ , where b is the harmonic oscillator 
size parameter, which has been taken from Ref. ^]. These widths go from = 3 MeV 
for f Si to = 20 MeV for ^He, as indicated in the upper panels of the just mentioned 
figures. 

b) IPSM-b: In the lower panels of the same figures are shown the results obtained when 
the recoil is convoluted with the Breit-Wigner distributions fl35|) for the strength func- 
tions of the fragmented deep hole states. The widths 7^-^^ have been estimated from 

12 



Refs. J26|, 123, l28|, l29|, l30|, IMI, l32|, l33|, [34, l35|, and in particular from [29, Figure 11] 



Table 1], with the results: 71S1/2 — ^ MeV in ]^C, 



= 14 MeV and 

7iP3/, = 3 MeV in fO \ and ^u,/, = 16 MeV and 71^3^, = 71^^^, = 5 MeV in f Si, 
both for protons and neutrons. One sees that, except for the ground states, the narrow 
peaks engendered by the recoil effect become now pretty wide bumps. 



We feel that the above rather rudimentary parameterization could be realistic enough for 
a qualitative discussion of the kinetic energy sum spectra. A more accurate model should 
be probably necessary for a full quantitative study and comparison with data. 
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FIG. 2: (Color online) Normalized energy spectra Snn{E)/T,^ for iHe, ^He, i^C, fO, and f Si hypernuclei 
for the full OMEP, obtained within the approaches IPSM-a (upper panel) and IPSM-b (lower panel). For 
the s-shell hypernuclei, only the IPSM-a approach has been used. 



^ The 3/2 peak is at 6.32 MeV, but small amounts of the P3/2 strength are also fragmented to the states 
of 9.93 MeV and 10.7 MeV 
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IV. GENERAL CONSIDERATIONS AND CONNECTION WITH DATA 



The normalized spectra shown in Figs. [T] and [2] have a very weak dependence on the 
dynamics involved in the NMWD process proper, and almost identical shapes would have 
been obtained if only the One-Pion Exchange Potential (OPEP) had been taken into account. 
To understand this fact one can appeal to the s-wave approximation, which assumes that only 
the relative matrix elements of the form (p, ISJT\V\0JJT), i.e., with the AA^ system in an s- 
state, significantly contribute to MjpPlLXSJT; j^). This has been examined quantitatively, 



for p- shell hypernuclei, in Ref. 



39| : see also the Refs. Furthermore, as we have 



discussed in Ref. [25|], those matrix elements depend only very weakly on the relative 
momentum p and, as such, they can be evaluated at the maximum value oi p = = V MA, 
which corresponds to P = or, according to Eq. fl20|) . to E = Aj^ = A. Thus the energy 
dependence of SnN{E)/TN remains exclusively in kinematical phase-space factors, and the 
position and width of the peaks will be unaffected by the dynamics of the decay process, 
which will influence, to some extent, only their relative heights. This is illustrated in the 
case of the np spectrum for in Fig. [31 showing that even this latter effect is very small. 
The comparison between the normalized spectra obtained with the full OMEP and with the 
7T + K exchange potential would give an almost perfect superposition. 

Needless to stress that the transition probabilities F^y do strongly depend on the hyper- 
nuclear transition potential, but this dependence is washed out in the ratios defining the 
normalized spectra. Conversely, for a given choice of transition potential both shell model 
approaches discussed here yield very similar results for the F^v. These points are illustrated 
in Table [U where one can also see that, as already anticipated in Section [Tll the effect of the 
recoil on these quantities is negligible. 

One can then summarize our findings by saying that in light systems (^He and \iie) 
the kinetic energy sum coincidence spectra SnN^E), normalized to the total decay rates Fjy, 
basically depend on energies associated with the three-body kinematics. The differences 
between Snp{E) and Snn{E) are mainly due to the differences in the proton and neutron 
separation energies and in the spectroscopic factors Fj^. For the remaining hypernuclei it 
is imperative, in addition, to take into account that most of the hole states are fragmented 
and consequently one has to consider the spreading of their strengths. 

The residual interaction among the valence particles and their coupling to the collective 
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TABLE I: Nonmesonic decay rates in units of = 2.50 x 10~^ eV and n/p branching ratios for 
j^^C, O and ^^Si computed with several transition potentials and using the IPSM-a and IPSM-b 
approaches. The values obtained in the IPSM-a framework but neglecting the recoil are shown 
within parentheses. 

A 

IPSM-a 

OMEP 0.249 (0.249) 0.956 (0.960) 0.260 (0.259) 

-K + K 0.244 (0.244) 0.755 (0.758) 0.323 (0.322) 
IPSM-b 

OMEP 0.246 0.947 0.260 

-K + K 0.241 0.748 0.322 

OPEP 0.142 1.004 0.141 

A ^ 

IPSM-a 
OMEP 
-K + K 
IPSM-b 
OMEP 
-K + K 



0.290 (0.290) 
0.287 (0.287) 

0.285 
0.282 



1.024 (1.027) 
0.811 (0.813) 

1.009 
0.799 



0.283 (0.282) 
0.354 (0.353) 

0.282 
0.353 



2804 
A 

IPSM-a 
OMEP 

-K + K 
IPSM-b 
OMEP 
-K + K 



0.348 (0.348) 
0.341 (0.341) 



0.345 
0.338 



1.163 (1.164) 
0.934 (0.935) 

1.123 
0.903 



0.299 (0.299) 
0.365 (0.365) 

0.307 
0.374 



15 




FIG. 3: (Color online) Normalized energy spectra Snp{E)/Tp for the decay of \^C, computed with the full 
OMEP and the OPEP in the framework of the IPSM-b. 

rotational and/or vibrational motions are not explicitly considered in the present work. But 
it is not expected that these would modify qualitatively the above scenario. It might be 
worth noting, nevertheless, that the pairing force would be capable of shifting some of the 
strength from the occupied levels to higher lying orbitals. For instance, in a part of the 
1^3/2 strength would be moved up into the lpi/2 orbital, while in |^Si the strength would 
be moved from the 1^5/2 level into the empty 2si/2 and 1^3/2 states. However, we believe 
that the coupling of the deep-hole states to other more complicated configurations through 
the residual interaction, which is treated here in a phenomenological way, is by far a more 
relevant effect for the physics discussed in the present paper. 

We would not like end this paper without making some comments on the relation between 
the formalism developed here and the experiments. To this end we shall follow closely the 



discussion in Ref. 



40|. The theoretical prediction for the number of nN pairs detected in 



coincidence with kinetic energy sum E within the interval dE can be written as 

(mnN{E) = CnN{E)SnN{E)dE, (37) 

where the factor CnN{E) depends on the experimental environment and includes all quanti- 
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ties and effects not considered in SnNiE), such as the number of produced hypernuclei, the 
detection efficiency and acceptance, etc. Assuming, for simphcity, that the experimental 
spectra have aheady been corrected for detection efficiency and acceptance and that the 
possible remaining energy dependence in this factor can be neglected, the predicted total 
number of detected events N„Ar can be related to F^v as follows: 

NnAT = J ^ ^-dE = CnN J SnN{E)dE = CnN'^N- (38) 

This allows us to rewrite (1371) in the form 

d^uNjE) SME) 

- ^^^^ 

What is measured in an experiment is the number of pairs AN^^{Ei) at a given energy 
Ei within a fixed energy bin AEnN, i-e., A'N'^^ (Ei) / AEnN- The total number of observed 
events is 

m 

N-^ = $:ANr^(E,), (40) 

4 = 1 

where m is the number of bins. The spectrum SnN{E) can be normalized to the experimental 
one by identifying in dSS]) with N^^. Thus, the quantity that we have to confront with 
measurements is 

A^ME) = K^AE,.N^4^. (41) 

tl 

For instance, to compare the experimental data given in 23!, Fig. 11 ] with our calculations 
shown in Figs. [1] and [21 the latter should be multiplied by the factors N^^^ ■ AEnp = 87 x 5 
MeV = 435 MeV, and N^^^ ■ AE„„ = 19 x 5 MeV = 95 MeV, respectively. 

Here we have, neither considered the resolution of the detector system, nor explicitly 
included the FSI. Also, we have ignored the three-body ANN decay contributions. As a 
consequence it is very reasonable that only a wide bump at about 140 MeV would appear in 
the experimental spectra for and heavier hypernuclei. The IPSM predicts quite similar 
spectra for the np and nn pairs. Moreover, from the present results one could venture to say 
that the neutron bump should lie at a smaller energy than the proton one. This agrees only 
marginally with the experiments performed so far, where important differences between the 
np and nn spectra have been observed. 
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V. CONCLUDING REMARKS 



In this paper, we have investigated, in the framework of the independent particle shell 
model (IPSM) , the effects of the recoil of the residual nucleus and of the spreading in strength 
of the deep-hole states on nonmesonic weak decay (NMWD) observables. We conclude that, 
while their effect is of minor importance for integrated observables like the decay rates 
and the n/p branching ratio, they play a crucial role in determining the shapes of the 
normalized kinetic energy sum coincidence spectra of nn and np pairs. For the spectra of 
s-shell hypernuclei, the recoil effect is the most important one. 

In summary, we believe that the IPSM is the appropriate lowest-order approximation for 
the theoretical calculation of the two-particle spectra in the NMWD when: 1) the recoil 
effect is included, and 2) the fragmentation of the strengths of the deep-hole states is taken 
into account. It is in comparison to this picture that one should appraise the effects of the 
FSI and of the two-nucleon-induced decay mode. 

The consequences of the two effects dealt with here on the one-particle kinetic energy 
spectra and on the opening angle distributions of np and nn pairs will be discussed elsewhere. 
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